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Abstract: In the present study, postnatal ontogenetic size and shape changes in the cranium of two lagomorph species, the plateau 
pika (Ochotona curzoniae) and woolly hare (Lepus oiostolus), were investigated by geometric morphometrics. The ontogenetic 
size and shape changes of their cranium exhibited different growth patterns in response to similar environmental pressures on the 
Qinghai-Tibetan Plateau. The overall size change in the cranium of the plateau pika was slower than that of the woolly hare. The 
percentage of ontogenetic shape variance explained by size in the woolly hare was greater than that in the plateau pika. The 
overall shape of the cranium was narrowed in both species, and morphological components in relation to neural maturity showed 
negative allometry, while those responsible for muscular development showed isometric or positive allometry. The most 
remarkable shape variations in the plateau pika were associated with food acquisition (temporalis development), though other 
remarkable shape variations in the incisive and palatal foramen in the ventral view were also observed. The most important shape 
change in the woolly hare was demonstrated by the elongation of the nasal bones, expansion of the supra-orbital process and 


shape variation of the neurocranium. 
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Leporidae, especially extant Lepus species, are am- 
ong the most capable mammalian runners. They can 
jump significant distances (up to 6 meters), move in 
rapid zigzag patterns, and locomote bipedally if required. 
The craniums of leporids are characterized by 
fenestration on the side of the rostrum that extends to the 
posterior region of the skull, an anatomical feature 
believed to reduce skull mass (Bramble, 1989). The 
posterior cranium in hares also permits considerable 
intracranial movement, which may function as a shock- 
absorbing mechanism to minimize jarring during high- 
speed running (Bramble, 1989). However, the extant 
Ochotonidae species from the genus Ochotona, which 
are closely related sister taxa of Leporidae, show less 
ability in running. When encountering predators, they 
usually warn other individuals through highly specialized 
calls and rapidly seek shelter. The shape of the incisive 
and palatal foreman in Ochotona have been used as 


diagnostic characteristics in taxonomy (Lissovsky et al, 
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2008; Smith & Xie, 2009). The body size evolution of 
Leporidae and Ochotonidae demonstrates markedly 
different patterns, with extant species of the former group 
showing non-normal distribution, and extant species of the 
latter group showing significant normal distribution (Ge et 
al, 2013). However, why or how size changes in each 
group correspond to the overall shape variance of the 
cranium during postnatal development is unclear. 

Of the known hares and pikas, the woolly hare (L. 
Oiostolus) and plateau pika (O. curzoniae) sympatrically 
inhabit a wide range of the Qinghai-Tibetan Plateau (QTP) 
and are both highly adapted to the cold and anoxic envir- 
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onment. Cranial growth patterns of these two species have 
not been studied, however, and it is unclear how different 
components of their cranium obtain their adult size and 
shape. 

In the present study, we aimed to explore the size 
and shape changes in the cranium of the plateau pika and 
woolly hare through postnatal ontogeny, and discern the 
possible causes for the divergence in the growth rates of 
different morphological components. Geometric morpho- 
metrics were used to quantify and visualize the morpho- 
logical variation in the dorsal, ventral and lateral cranium. 


MATERIAL AND METHODS 


Materials 

A total of 84 cranium of L. oiostolus and 56 speci- 
mens of O. curzoniae were investigated. All cranium 
specimens were intact on at least one side of the cranium. 
The specimens represent the two species at different ages, 
and were mainly collected from the Qinghai-Tibet 
Plateau of China. All materials were preserved in the 
National Zoological Museum of China, Institute of Zool- 
ogy, Chinese Academy of Sciences (IOZCAS) and the 
Northwest Institute of Plateau Biology, Chinese Acad- 
emy of Sciences (NIPBCAS). Previous studies demonst- 
rated insignificant sexual dimorphism in Ochotona and 
Lepus, which enabled the combination of males and 
females in all our analyses (Lu, 2003; Smith & Weston, 
1990). Supplementary details on the specimens studied 
are given as Supplementary Table 1 (Supporting inform- 
ation of http://www.zoores.ac.cn/). 


Photography 

We viewed the dorsal, ventral and lateral cranium in 
this study. Photographs were taken using a Canon Power- 
Shot SSIS (Japan) with a macro-focusing lens. To mini- 
mize measurement errors, all photos were taken by the 
corresponding author. For dorsal views, the camera was 
focused on the frontal region; for ventral views, the 
camera was focused on the tooth rows; and for lateral 
views, the camera was focused on the zygomatic arch. 
Images were standardized for skull position and camera 
lens plane as well as for distance between camera lens 
and the sample. All digital photos were taken with a scale 
bar parallel to the longitudinal axis of the cranium, which 
was set to record the size differences among specimens. 


Data acquisition and partition 
In total, 129 landmarks and semi-landmarks on the 
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right half of the skull in O. curzoniae were digitized for 
analysis, including 39 for the dorsal cranium, 55 for the 
ventral cranium, and 35 for the lateral cranium (Figure 
1A1, A2 and A3). A total of 168 landmarks and semi- 
landmarks on the right half of the skull in L. oiostolus 
were digitized for analysis, which included 44 for the 
dorsal cranium, 62 for the ventral cranium, and 62 for the 
lateral cranium (Figures 1B1, B2 and B3). The locations 
of these landmarks and semi-landmarks mainly followed 
those of Ge et al (2012); detailed information is given in 
Table 1. Terminologies used in this study mainly 
followed Wible (2007). The landmarks and semi- 
landmarks were defined and digitized in tpsDig 2.0 
(Rohlf, 2001) with MakeFan6 (Sheets, 2003), which 
offered identifiable and homologous locations on the 
skull. The scale factor was also generated in the data 
matrix with tpsDig, which set 1 cm of the scale bar 
recorded in the original images as 1. The raw datasets of 
each of the three views were examined in tpsSmall to test 
whether any specimens deviated largely from the whole 
dataset. Semi-landmarks were permitted to slide along 
the tangents to the curve to minimize bending energy 
between each specimen. The skulls of mammals can be 
divided into pre-orbital, inter-orbital and post-orbital 
regions (Lin & Shiraishi, 1992; Tanaka, 1942), which 
were extracted from the ventral view in the present study 
to respectively compare their ontogenetic size changes. 


Morphometric analysis 

Morphometric analysis of ontogenetic size and shape 
change in O. curzoniae and L. oiostolus was conducted 
using a generalized orthogonal least-squares Procrustes 
superimposition and relative warps analysis. Statistical 
analysis was based on centroid size, that is, the square 
root of the sum of squared distances of measured land- 
marks. Visualization of shape change was performed by 
thin-plate spline procedures. These analyses were condu- 
cted in the tps series of programs (Rohlf, 2001) and 
SPSS statistics 17.0 (SPSS inc., Chicago, IL, USA). 

In the present study, centroid size was calculated 
from the ventral view of the whole cranium, and the pre- 
orbital, inter-orbital, and post-orbital regions in each of 
these two species. Growth rate (inferred from size chan- 
ge) of these regions was compared respectively. Centroid 
sizes of each of the three regions were regressed against 
the centroid size of the whole cranium to investigate the 
relative growth rate of each morphological partition. The 
null hypothesis was the size change of these small 
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partitions parallel to the size change of the whole cran- region in the plateau pika grew slightly slower than that in the 
ium. The coordinate data of each view were regressed woolly hare, while the post-orbital region of the woolly hare 
against centroid size in tpsRegr (Rohlf, 2001) to showed a no-linear growth pattern (P=0.086) (Figure 2). 
visualize the ontogenetic trajectories in different regions Regressing size change on shape change demonstr- 
of the cranium by the direction and length of vectors. ated prominent allometry in the three views of the woolly 
The Generalized Goodall F-test of significance (Goodall, hare (dorsal: (—43.37, df=84, 6384, P<0.01, size 
1991), with dfl and df2 degrees of freedom, was cond- explained 36.03% of total shape variation; ventral: 
ucted. Shape changes from the smallest specimen to the FH 26.2383, df=120, 9840, P<0.01, size explained 24.13% 
largest specimen were illustrated by deformation grids. of total shape variation; lateral: /=26.87, df=120, 8400, 
P<0.01, size explained 27.24% of total shape variation), 
RESULTS and in the dorsal and ventral views of the plateau pika 


(dorsal: /=10.82, df=74, 3996, P<0.01, size explained 
Analyses in tpsSmall indicated that no specimen 16.64% of total shape variation; ventral: 15.38, df=104, 
deviated from other specimens in each view for the two 5512, P<0.01, size explained 22.43% of total shape 


species (r ranged from 0.9999 to 1 for L. oiostolus and O. variation), but a lower value was observed in the lateral 
curzoniae in tpsSmall), which allowed us to use infant view of the plateau pika (0.25, df=66, 3234, P<0.01, 
specimens in the present study. Geometric morphometric size explained 0.5% of total shape variation). These 
analyses demonstrated allometric growth in different cran- results indicated that the percentage of total ontogenetic 
ium regions of the plateau pika and woolly hare. The crani- shape variance explained by size in the woolly hare was 
um growth rates inferred from size change in the pro-orbital, higher than that in the plateau pika. 

inter-orbital, and post-orbital regions were higher in the Plotting shape change vectors on consensus shape 
woolly hare, though the growth of the pro-orbital region demonstrated that shape variation in the two species 
was nearly parallel in the two species. The inter-orbital mainly occurred in the very early life stage (Figure 3). 





Figure 1 Landmark and semi-landmark locations 


A1-A3: Dorsal, ventral and lateral views of O. curzoniae; B1-B3: Dorsal, ventral and lateral views of L. oiostolus. 
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Figure 2 Size variation (centroid size) in the pre-orbital, inter-orbital and post-orbital regions of the cranium 


Table 1 Definition of landmarks and semi-landmarks 


O. curzoniae 


L. oiostolus 





Dorsal cranium 

1. Anterior tip of nasal bone 

2-10. Semi-landmarks along lateral margin of premaxilla and maxillary 

11. Anterior tip of zygomatic in the lateral edge 

12-17. Semi-landmarks along lateral edge of zygomatic arch and squamosal 

18. Most posterior end of zygomatic arch 

19. Meeting point between nasal bone and frontal along longitudinal axial of 
cranium 

20. Posterior tip of nasal bone 

21. Most anterior point of orbital 

22. Mid point between 21 and 23 

23. Most internal point of orbital 

24. Mid point between 23 and 25 

25. Most posterior point of orbital 

26-28. Semi-landmarks along internal edge of zygomatic arch 

29. Posterior point of the zygomatic process of squamosal body 

30. Posterior point at the joining of squamosal body to the zygomatic process of 
squamosal 

31-35. Semi-landmarks along posterior edge of cranium 

36. Meeting point of frontal and parietal along longitudinal axial of cranium 

37. Most anterior meeting point of frontal and parietal 

38. Meeting point of parietal and supraoccipital along longitudinal axial of cranium 

39. Most posterior point of the supraoccipital along longitudinal axial of cranium 

Ventral cranium 

1.Meeting point of incisors 

2. Most posterior end of incisors along lateral edge 

3-10. Semi-landmarks along lateral edge of premaxilla and maxillary 

11. Most anterior edge of incisive and palatal foramen 

12-21. Semi-landmarks along incisive and palatal foramen 

22. Most posterior end of the incisive and palatal foramen 

23. Posterior end of the incisive and palatal foramen along longitudinal axial of the 
cranium 

24. Anterior end of tooth row 

25. Posterior end of tooth row 

26. Most anterior end of masseteric spine 

27. Apex of masseteric spine 

28. Most posterior end of masseteric spine 

29-31. Semi-landmarks along lateral edge of zygomatic arch and the zygomatic 

process of squamosal 

32. Most posterior end of the zygomatic process of squamosal 

33. Most anterior end of orbital 

34. Most internal end of orbital 

35. Most posterior end of orbital 

36-37. Semi-landmarks along lateral edge of orbital 

38. Posterior end of palatine along longitudinal axis 

39-40. Semi-landmarks along entopterygoid crest 

41. Posterior point of zygomatic process of squamosal body 

42. Contact point of glenoid fossa with alisphenoid 

43. Mid point of basioccipital 

44. Anterior point of tympanic bulla 

45. Posterior point of tympanic bulla 

46-48. Semi-landmarks along lateral edge of tympanic bulla 

49-52. Semi-landmarks along internal edge of tympanic bulla 
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Dorsal cranium 
1. Anterior tip of nasal bone 
2-12. Semi-landmarks along lateral margin of premaxilla and 
maxillary 
13. Anterior tip of zygomatic in the lateral edge 
14-18. Semi-landmarks along lateral edge of zygomatic 
19. Meeting point between nasal bone and frontal along longitudinal 
axial of cranium 
20. Mid point between landmark 19 and 21 
21. Posterior tip of suture between nasal bone and frontal 
22-26. Semi-landmarks along interior edge of zygomatic 
27. Anterior meeting point between zygomatic and frontal 
28. Mid point between 27 and 29 
29. Most posterior beginning of anterior supraorbital process 
30. Anterior tip of anterior supraorbital process 
31-32 Semi-landmarks between landmark 30 and 33 
33. Posterior tip of posterior supraorbital process 
34. Most anterior beginning of posterior supraorbital process 
35-36. Semi-landmarks along posterior contour of orbital 
37. Meeting point of frontal and parietal along longitudinal axial of 
cranium 
38. Meeting point of parietal and supraoccipital along longitudinal 
axial of cranium 
39. Most posterior point of the supraoccipital along longitudinal 
axial of cranium 
40. Most posterior point of parietal 
41. Most anterior meeting point of parietal and supraoccipital 
42. Most lateral point of parietal 
43. Semi-landmarks along lateral margin of squamosal 
44. Most anterior point of parietal meeting with frontal and 
squamosal 
Ventral cranium 
1.Meeting point of incisors 
2-11.Semilandmarks along lateral edge of premaxilla and maxillary 
12.Most anterior point of zygomatic 
13-20. Semi-landmarks along lateral edge of zygomatic 
21-30. Semi-landmarks along internal edge of premaxilla and 
maxillary 
31. Anterior end of the incisive and palatal foramen along 
longitudinal axial of cranium 
32. Posterior end of the incisive and palatal foramen 
33. Anterior end of the cheek tooth row 
34. Posterior end of the cheek tooth row 
35. Posterior end of palatal bridge along longitudinal axial of 
cranium 
36. Anterior end of palatal bridge 
37-39. Semi-landmarks along lateral margin of entopterygoid crest 
40. Most anterior end of orbital 
41-45. Semi-landmarks along internal edge of zygomatic 
46. Anterior end of basisphenoid along longitudinal axial of cranium 
47. Most anterior point of basisphenoid meeting with basisphenoid 
48. Lateral end of basisphenoid meeting with basisphenoid 
49.Most anterior end of foreman magnum 
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O. curzoniae 


(Continued) 


L. oiostolus 





53. Most anterior point of foramen magnum 
54. Most posterior end of foramen magnum 
55. Most lateral point of foramen magnum 





Lateral cranium 

1.Center of alveolar ridge over maxillary incisor 

2-3. Semi-landmarks along lower edge of premaxilla and maxillary 
4. Most upper point of maxillary 

5. Most anterior end of tooth row 

6. Most posterior end of tooth row 

7. Most anterior end of check tooth 

8. Most anterior end of nasal bone 

9-18. Semi-landmarks along upper ridge of the cranium 

19. Posterior end of the cranium 

20.Most upper point of tympanic bulla opening 

21.Most lower point of tympanic bulla opening 

22. Most basal point of basioccipital 

23. Most lower end of tympanic bulla 

24. Most internal end of tympanic bulla 

25.Most lateral end of zygomatic arch 

26-29. Semi-landmarks along lateral edge of zygomatic arch 
30. Most posterior end of orbital 

31-35. Semi-landmarks along upper edge of zygomatic arch 
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Details are also shown in Figure 1. 
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50.Most posterior end of foreman magnum 
51. Most lateral end of foreman magnum 

52. Opening of external carotid artery 

53—62. Semi-landmarks along tympanic bulla 





Lateral cranium 

1.Center of alveolar ridge over maxillary incisor 

2-11. Semi-landmarks along lower edge of premaxilla and 
maxillary 

12. Anterior end of cheek tooth raw 

13. Posterior end of cheek tooth raw 

14. Most anterior point of zygomatic 

15-24. Semi-landmarks along lower edge of zygomatic 

25. Most posterior end of zygomatic 

26. Most posterior end of orbital 

27-40. Semi-landmarks along orbital 

41. Anterior end of supraorbital process 

42. Posterior end of supraorbital process 

43. Anterior end of incisive 

44-47. Semi-landmarks along upper edge of premaxilla and 
maxillary 

48. Anterior end of nasal bone 

49-51. Semi-landmarks along upper edge of nasal bone 

52-53. Semi-landmarks along posterior cranium 

54. Most superior end of occipital 

55. Posterior end of occipital 

56. Most lower end of paroccipital process 

57 Most lower point of tympanic bulla opening 

58 Most upper point of tympanic bulla 

59 Most anterior point of tympanic bulla 

60. Most lower point of tympanic bulla 

61. Meeting point of frontal, parietal and sphenoid bones 

62. Meeting point of supraoccipital, parietal and sphenoid bones 
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Figure 3 Growth trajectories in different views of the cranium 
Al-Aé6: O. curzoniae, Al—A2, dorsal view; A3—A4, ventral view; A5—A6, lateral view. B1-B6: L. oiostolus, B1—B2, dorsal view; B3—B4, ventral view; 


B5-B6, lateral view. 1, 3, 5 of AB show variation in the early life stage; 2, 4, 6 show variation at near adulthood. 
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The overall shape of the plateau pika was narrowed, and 
the anterior ridge of the zygomatic arch was expanded 
slightly. Incisive and palatal foreman showed prominent 
shape variation from infants to adults. The lateral view of 
the cranium showed less shape variation compared with 
the other two views. The woolly hare exhibited a wider 
range of shape variation than that of the plateau pika, and 
the general shape of its cranium was also narrowed with 
a remarkable shape change in the nasal bone region and 
contraction of the tympanic bulla. The super-orbital 
process displayed prominent enlargement in both the 
anterior and posterior branches. The relative size of the 
neurocranium was reduced, the relative size of the tymp- 
anic bulla was greatly reduced and the incisive and palatal 
foramen showed less variation. For plateau pika, regre- 
ssing centroid size on shape variables also indicated that 
small morphological components in the pre-orbital region 
(in all three views) was dominated by positive allometry, 
while those in the inter-orbital and post-orbital regions 
(dorsal and ventral views) showed prominent negative 
allometry. In the woolly hare, however, morphological 
eames Positive allometry 
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components in the pro-orbital and inter-orbital regions 
were mainly dominated by isometry and positive allom- 
etry, while a prominent negative allometry was observed 
in the post-orbital region. A detailed comparison of the 
shape variation in small morphological components 
between the two species is given in Figure 4. 


DISCUSSION 


In light of the highly specialized skull morphology 
and rapid skull development in the plateau pika and 
woolly hare, as well as their short period of offspring 
dependence, the special growth patterns of their cran- 
iums are likely the result of combined pressure for food 
acquisition and other strategies for survival. Natural 
selection may favor the rapid development of adult 
morphology to improve food acquisition, locomotion pe- 
rformance, and predator avoidance. The occurrence of 
different growth rates in specific morphological comp- 
onents suggests selective pressures on different parts of 
the cranium during postnatal growth. 
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Figure 4 Shape variation from the smallest to largest specimens 
Al—A6: O. curzoniae, A1—A2, dorsal view; A3—A4,ventral view; A5—A6, lateral view. B1—B6: L. oiostolus, B1—B2, dorsal view; B3—-B4, ventral view; 


B5-B6, lateral view. Arrows indicate prominent shape variation in incisive and palatal foreman of O. curzoniae and supra-orbital process of L. oiostolus. 
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